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Abstract — Broadcasting systems have to deal with channel 
diversity in order to offer the best rate to the users. Hierarchical 
modulation is a practical solution to provide several rates in 
function of the channel quality. Unfortunately the performance 
evaluation of such modulations requires time consuming simu- 
lations. We propose in this paper a novel approach based on 
the channel capacity to avoid these simulations. The method 
allows to study the performance in terms of spectrum efficiency of 
hierarchical and also classical modulations combined with error 
correcting codes. Our method will be applied to the DVB-SH 
standard which considers hierarchical modulation as an optional 
feature. 

I. Introduction 

In most broadcast applications, the receivers do not expe- 
rience the same signal-to-noise ratio (SNR). For instance, in 
satellite communications the channel quality decreases with 
the presence of clouds in Ku or Ka band or with shadowing 
effects of the environment in lower bands. 

The first solution for broadcasting is to design the system for 
the worst-case reception. However this solution does not take 
into account the variability of channel qualities. This holds 
a loss of spectrum efficiency for users with good reception. 
Then two other schemes have been proposed in [T] and 
to improve the first one: time division multiplexing with 
variable coding and modulation, and superposition coding. 
Time division multiplexing consists in using a first couple 
modulation/coding rate during a fraction of time, and then 
using an other modulation/coding rate for the remaining time. 
All the population can receive the first part of the signal called 
the High Priority (HP) signal and only the receivers in good 
conditions receive the second one called the Low Priority (LP) 
signal. 

Unlike time sharing, superposition coding sends information 
for all the receivers all the time and it was proved to achieve 
better rates than time sharing for the continuous Gaussian 
channel fl2). In superposition coding the available energy is 
shared to several service flows which are sent simultaneously 
and in the same band. Hierarchical modulation is a practical 
implementation of superposition coding. Figure [T] presents the 
principle of the hierarchical modulation using a non-uniform 
16-QAM. The idea is to merge two different streams at the 
modulation step. The HP stream is used to select the quadrant, 



and the LP stream selects the position inside the quadrant. In 
good conditions receivers can decode both streams, unlike bad 
receivers which only locate the quadrant and then decode the 
HP stream as a QPSK constellation. 
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Fig. 1: Hierarchical Modulation using a non-uniform 16-QAM 

In practical contexts, both schemes can be used com- 
bined with applications generating hierarchical flows, like 
H.264/SVC video 0, 0. 

Satellite communications standards like DVB-S2 [5| and 
DVB-SH [6] consider hierarchical modulation as an optional 
feature. To our knowledge, the only way to estimate the 
performance of real code using hierarchical modulation is to 
run simulations. We propose here a new approach to evaluate 
the performance of such modulation. The method is based on 
the channel capacity and relies on the fact that the real code 
at coding rate R works as a theoretical ideal code at coding 
rate R. Our approach is applied to DVB-SH. 

The paper presents our work as follows: 

• In Section|II]we compute the capacity for any hierarchical 
modulation. 

• We propose a method using the capacity to evaluate 
the performance of hierarchical modulations in terms of 
spectral efficiency and required E s /Nq for a targeted 
Bit Error Rate (BER) or Packet Error Rate (PER) in 
Section Hill 



Section IV concludes the paper by summarizing the results 
and presents the future work. 



II. Capacity of the Hierarchical Modulation 

In this section we introduce the channel capacity and 
compute it for any hierarchical modulation. Some results are 
given at the end of this section using the hierarchical 16-QAM. 

A. Computation of the capacity 

A channel can be considered as a system consisting of 
an input alphabet, an output alphabet and a transition matrix 
p(y\x). We define two random variables X, Y representing the 
input and output alphabets respectively. The mutual informa- 
tion between X and Y, noted I(X; Y), measures the amount 
of information conveyed by Y about X. For two discrete 
random variables X and Y, ([T]i gives the expression of the 
mutual information. It can be extended to continuous random 
variables 0. 
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Using this notion, the channel capacity is then given by 

C = maxI[X;Y), (2) 

p(x) 

where the maximum is computed over all possible input 
distributions. 

Here we consider the memoryless discrete input and con- 
tinuous output Gaussian channel. The discrete inputs Xj are 
obtained using a modulation and belong to a set of discrete 
points x C K 2 of size \x\ = M = 2 m called the constellation. 
Thus, each symbol of the constellation carries m bits. An 
explicit formula (|3]l for the capacity in this particular case 
is given in [8, Chapter 3]. 



C = log 2 (M)- 
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We are now interested to evaluate the capacity of a stream 
using a subset of all the bits. Reference [9| presents the 
case where each stream uses one bit. The idea is to modify 
the input random variable X by the input bits used in each 
stream. We define bj as the value of the i th bit of the label 
of any constellation point x. Suppose the stream uses k bits 
among m in the positions i±, For any integer i, let 

l n {i) denote the n th bit in the binary representation of i such 
as i — J2n=i 'n(i)2 n_1 - Then using M with the new input 
variable and continuous output, the capacity of the stream is 
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where p(y) = A J2 i=0 P(y\K = h(i), -M k = **(«)). 



Let L n (x) denote the n th bit of the label of x. We introduce 
Xi the subset of x defined as follows 



(5) 



The set Xi depends on i and the positions of the bits involved 
in the stream. Figure [2] shows an example of subsets for a 16- 
QAM with a particular mapping, where the stream uses bits 
in position 1 and 2 (in that case k = 2, i\ = 1 and i 2 = 2). 
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Fig. 2: Examples of xo (vertical lines) and X3 (horizontal lines) 

Then the conditional probability density function of y in 
Q can be written 



p (y\K = h(i), -,b ik = l k (i)) = p(yl x M x l x e x%) 



where \x%\ = 2 m k for all i. Moreover, the transition 
distribution p(y|x) for a Gaussian channel is 
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Using |6]) and |7]i in Q, we finally obtain the capacity 
for one stream The capacity is an increasing function of 
E s /N and its value is less or equal to k the number of bits 
used in the stream. The positivity of the capacity results of its 
mathematical definition Q. 
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B. Case of non-uniform hierarchical modulation capacities 

We begin with few definitions before applying to the 
hierarchical 16-QAM considered in DVB-SH. 

Hierarchical modulations merge several streams in a same 
symbol. They often use non-uniform constellation. Non- 
uniform constellations are opposed to uniform constellations 
where the symbols are uniformly distributed. The constellation 
parameter is defined to describe non-uniform constellations. 
Figure [T] illustrates a non-uniform 16-QAM. The constellation 
parameter a is defined by a = d^/di where 2dh is the 
minimum distance between two constellation points carrying 
different HP bits and 2d; is the minimum distance between 
any constellation points. Thus a verifies a > 1 where a = 1 
corresponds to the uniform 16-QAM. The DVB-SH standard 
[10 1 recommends two values for a: 2 and 4. 

Back to the hierarchical modulation capacity, we suppose 
the HP stream uses the bits in position 1 and 2. The LP stream 
involves the remaining bits, i.e., bits 3 and 4. Figures [3] and [4] 
present the 16-QAM capacity for the two values of a defined 
in the DVB-SH guidelines ifTOl - 

When a grows, the constellation points in each quadrant 
become farer to the I and Q axes. Thus it is easier to decode the 
good quadrant and the capacity for the HP stream increases. 
However in the same quadrant, the points gets closer and the 
capacity for the LP stream decreases as it is harder to decode 
which symbol was sent. Concerning the HP stream capacity, 
it increases with a but nevertheless it is limited by the QPSK 
capacity as the same energy is used for 16 points instead of 
4. 

To conclude, we can observe that the global hierarchical 
modulation capacity (i.e. HP+LP) is always less than the 16- 
QAM capacity for any value of a > 1. This result has been 
proved in [9] where the different streams use one bit. We can 
also note that the best total capacity is achieved when the 
constellation is close to the uniform constellation, but on the 
other hand performances for the HP stream are decreased. 




Fig. 3: 16-QAM Hierarchical Modulation Capacity - a = 2 

III. Performance Evaluation for Hierarchical 
Modulations 

We present in this part a method to get easily the spec- 
trum efficiency and the required E s /Nq for a given targeted 
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Fig. 4: 16-QAM Hierarchical Modulation Capacity - a = 4 



BER/PER for any modulation and coding rate without comput- 
ing extensive simulations. In [11] a fast coding/decoding per- 
formance evaluation method based on the mutual information 
computation is proposed and applied to time varying channel 
for QPSK modulation with very good prediction precision. 
The method developped hereinafter based on channel capacity 
computation makes it possible to predict the performances of a 
coding scheme combined with any modulation and especially 
any hierarchical modulation. 

A. Principle 

Before applying the approach to a real code, we begin with 
an example using theoretical ideal codes achieving the channel 
capacity. The normalized capacity for a modulation is defined 
by C mod = ^C mod where C mod is the modulation's capacity 
and m denotes the number of bits per symbol. The C mod 
function has some properties: 

• It is an increasing function of E s /Nq. 

• It verifies < C mod < 1. 

• It is equivalent to the notion of coding rate when using 
ideal codes. 

Given a modulation and a coding rate, if we want to know 
at which E s /Nq the code is able to decode, we just need to 
inverse the normalized capacity function for that coding rate. 
Figure [5] illustrates this example using a QPSK. 
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Fig. 5: Decoding threshold for a coding rate of 0.8 



Let us now consider an actual waveform based on a real 
code. We can consider that the real code at a coding rate R 
works as an ideal code with rate R where R> R. Determining 
R requires a performance curve (BER or PER vs E s /Nq) 
for one actual reference modulation and coding scheme with 
code rate R. The performance curve allows to estimate the 
behaviour of the code. We illustrate our method with an 
example. We would like to study the hierarchical modulation 
using an non-uniform 16-QAM {a = 2) in the DVB-SH 
standard at a target BER of 10~ 5 . We use the 2/9 and 1/5-turbo 
codes for the HP and LP streams respectively. To determine 
R and the decoding thresholds, the method works as follow: 

1) Use the performance curve of the reference modulation 
with rate R to get the operating point (E s /No) re f 
corresponding to the desired performance. In the DVB- 
SH guidelines [10, Table 7.5], we read: 

For the coding rate 2/9, BER QPSK (-3.4dB) = 1(T 5 
For the coding rate 1/5, BERq PSK (-3.9dB) = 1CT 5 

2) Compute the normalized capacity for the reference mod- 
ulation which corresponds to R: 

For the HP stream, R = <7q PSK (-3.4dB) « 0.27 
For the LP stream, R = C QPS k (-3.9dB) « 0.25 

3) For the studied modulation, compute E s /Nq such as the 
normalized capacity at this SNR equals R: 

(E s /N ) Hp = C^ a=2 (A = 0.27) = -2.7dB 

(E s /N ) Lp = C^p a=2 (R = 0.25) = 6.2dB 

The DVB-SH guidelines provide the decoding thresholds 
for the HP and LP streams ifTUl Figure 7.40]. Less the 
0.3dB due to the pilots, we read -2.6dB and 6.2dB for 
the HP and LP streams respectively. 

4) Finally the points ((E S /N ) HP , Rhp X m) and 
((E s /Nq) lp , Rlp x tti) are plotted on the spectrum 
efficiency curve (e.g., Figure [7] and |HJ. 

These steps are repeated for all the coding rates. Our 
method makes two assumptions. First of all, it approximates 
the information rate by R x m. This approximation is justified 
by the fact that the targeted performance (BER = 10~ 5 ) is 
very small and thus only hardly impact the useful information 
rate. 

The second assumption is to suppose as in ifTTl that the 
performance of the decoding only depends on the normalized 
capacity and not on the modulation as for ideal codes. 
To validate this second assumption, Figure [6] presents the 
normalized capacity in function of the coding rate for various 
modulations using the data of DVB-SH OH Table 7.5]. 

B. Application to DVB-SH 

In our study, we use the guidelines of DVB-SH 1101 Table 
7.5] to get the reference curves and the operating points at a 
target BER of 10 -5 . These data correspond to a static receiver. 
The guidelines also provide all the reference numerical results 
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Fig. 6: DVB-SH: normalized capacities 

for the hierarchical modulation IfTUl Figure 7.40] where 0.3dB 
need to be removed due to the pilots. Thus it allows to evaluate 
the efficiency of our method. 

The method described earlier is now applied to plot the 
spectrum efficiency curves as a function of required E s /Nq 
(target BER = 10~ 5 ). Figures [7] and [8] present the results for 
a = 2 and a = 4 respectively. The reference results from the 
guidelines correspond to the standard curves. The results show 
a good precision and it does not require computing extensive 
simulations. Moreover, our work consolidates the fact that the 
capacity is a good metric to evaluate performance. 




Fig. 7: DVB-SH spectrum efficiency, a = 2, BER = 10~ 5 
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Fig. 8: DVB-SH spectrum efficiency, a = 4, BER = 10~ 5 



C. DVB-SH: Additional Results 



We apply hereafter the method to compute the required 
E s /Nq in function of a or the coding rate. 

Figures [9] and 10 present how the required E s /N varies 
with a. For the LP stream, the required E s /N is an increasing 
function of a, unlike the one for the HP stream who decreases. 
In fact the required E s /Nq for the HP stream tends to the 
required SNR of the QPSK modulation. These results are 
obvious when we look the modification of the constellation 
with a. When a increases, the points in one quadrant become 
closer and the constellation is similar to a QPSK. It explains 
why the LP stream is harder to decode and requires a better 
SNR. 



Fig. 1 1 : Required E s /N a function of the coding rate, a = 2 



Fig. 9: Required E s /N function of a, HP Stream, BER 
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IV. Conclusion 

In this paper, we introduce a general method allowing to 
analyse the performance of hierarchical modulations. This 
method relies on the channel capacity which has been com- 
puted for any kind of constellation. It has been applied to 
DVB-SH in order to obtain the spectrum efficiency. The 
comparisons with reference simulation results show the good 
reliability of the method. 

Future work will improve the comparison between hier- 
archical modulation and time sharing. We also plan to take 
into account the mobility in the context of DVB-SH. Finally 
we will apply our method to DVB-S2 and other standards 
including hierarchical modulation to complete the study. 
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Fig. 10: Required E s /N function of a, LP Stream, BER 
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The last point concerns the variation of the required E s /N 
with the coding rate. It is presented on Figure 11 for a = 
2. Obviously, if there is less redundancy, the SNR has to be 
higher to decode. 



All these results are helpful during the dimensioning of a 
broadcast system. If the operator knows the SNR distribution 
of the receivers, it is possible to optimize the choice of the 
constellation parameter or the coding rate for each stream. 
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